Real-time vector field tracking with a cold-atom magnetometer 
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We demonstrate a fast three-axis optical magnetometer using cold, optically-trapped 87 Rb gas as a sensor. 
By near-resonant Faraday rotation we record the free-induction decay following optical pumping to obtain 
the three field components and one gradient component. A single measurement achieves shot-noise limited 
sub-nT sensitivity in 1 ms, with transverse spatial resolution of ^20 [J.m. We make a detailed analysis of the 
shot-noise-limitcd sensitivity. 
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Control of magnetic fields is critical to many ap- 
plications, for example high-sensitivity instruments 1 
and atomic physics experiments^. Hot-atom optical 
magnetometers^ offer sub-fT sensitivity at the cm scale^ 
and sub-pT at the mm scaled Cold atom magnetome- 
ters have d emo nstrated 10 pT gradient sensitivity at the 
10 \xm scald^ and 100 pT gradient sensitivity at the 
3 \im scaled A single-ion narrow-band (2 x 10~ 3 Hz) 
magnetometer^ showed 10 pT sensitivities. Many of 
these systems are designed for sensing of one field or 
gradient component. In contrast, precise control of 
fields requires simultaneous knowledge of all field compo- 
nents and possibly gradients as well. Recent work with 
hot-atom magnetometers has demonstrated three- axis 
sensing and control with nT precision at cm scales jSEol 
while a three-axis modulated cold-atom magn etometer 
has shown nT sensitivity at 300 (i.m scalesP^D2] 

Here we demonstrate sub-nT sensitivity at 50 (i.m 
length-scale in a cold-atom magnetometer employing 
near-resonant Faraday rotation probing^. The instru- 
ment gives three-axis field information plus one gradient 
component, obtained by free-induction decay (FID) af- 
ter optical pumping. A full measurement requires 1 ms, 
and can be repeated with zero dead-time, allowing high- 
bandwidth recording and control of the vector field. No 
external magnetic fields need be applied, so the technique 
can be used for real-time monitoring during field-sensitive 
processes. Our implementation is shot-noise limited and 
we give explicit expressions for its shot-noise-limited sen- 
sitivity. 

The experiment is shown schematically in Figure [lja) . 
An ensemble of atoms is held in an elongated optical 
trap, and subject to an unknown field B. The atoms are 
first optically pumped along the z direction, so that the 
collective atomic spin F achieves the value (0,0, F z (0)), 
and then allowed to precess around B. Off-resonance 
Faraday-rotation probing measures the rotation angle 
9 = GF Z , where G is a coupling constant, known from in- 
dependent measurements.^ We observe the FID signal 
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where 7 = hbSf^ is the gyromagnetic ratio, /is is 



the Bohr magneton, gp is the Lande factor, and h is 
Planck's constant. The transverse relaxation time T2 = 
l/(I>y| .By |) is due to the field-parallel gradient compo- 
nent £?'| = d\B\/d z , and a Lorentzian distribution (full- 
width at half-maximum L) of atoms along z, the trap 
axisP^ The process is then repeated with the spins ini- 
tially polarized (0, F y (0), 0), to give 
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Fitting the composite signal from the two FID mea- 
surements gives the three B components up to a global 
sign and T2. The ambiguity can be lifted by applying 
a known field, if necessary. Representative traces are 
shown in Figure [TJ c) . Re lative to other vector magne- 
tometry techniques j 9 l 10 l 12 l this method is simple both in 
procedure and in interpretation and requires no applied 
B-fields, making it attractive for work with field-sensitive 
systemsPES 

To derive Equations |l]) and ^ we note that the mi- 
croscopic spin operators evolve as P(t) — R(zi,t)P(0), 
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is the generator of rotations about B and B = B/|B|. 

Possible decoherence mechanisms include atomic mo- 
tions and collisions, tensorial light shifts due to the probe 
light, and decoherence due to gradient of the field. In 
our experiment the effect of tensorial term of probe is 
negligible, since we are far detuned from D2 transition 
line and we are using few photons for detection. For the 
time-scales involved in this experiment, decoherence due 
to collisions is negligible, whereas dephasing, i.e. differ- 
ential precession due to field inhomogeneity, typically is 
not. In the language of magnetic resonance, we expect 
longitudinal relaxation to be much slower than transverse 
relaxation due to field inhomogeneity. 
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Expanding the field as B(z) w B + (Bj| + B' ± )z + 
0(z 2 ), where Bji is parallel to Bo and B^ is perpendic- 
ular. We note that a change in the magnitude of B has 
an accumulating effect on the spin precession, i.e., the 
change in f grows with t. In contrast, a change in the 
direction of B has a fixed effect: From the perspective 
of the measurement, a rotation of B is equivalent to a 
rotation of both the initial state and the measured com- 
ponent F z . For small gradients d z B -C B/^ atoms , where 
'atoms is the length of the cloud, we can ignore B^ . This 
approximation, along with the fact that A B +2 = —A 
allows us to write 
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R(z,t) «I+ A Bo s'mw(z)t + A% [1 - cosu(z)t] (4) 

where w(z) = 7f|B + zB'^\. 

In our trap, we observe an atomic density p(z) well 
approximated by a Lorenzian p(z) — L/ir(z 2 + L 2 ) where 
L s» 48 pm is the full- width half-maximum extent of the 
ensemble. The collective spin F = J2i P then evolves as 
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In the first term I + A 2 B describes a projector onto the 
direction of Bo- This is the steady-state polarization. 
The second line describes a decaying oscillation of the 
transverse components, i.e., those perpendicular to Bo- 

Our experimental apparatus has been described in de- 
tail elsewhere.^ Briefly, we work with ensemble of up to 
2 x 10 5 laser cooled 87 Rb atoms in the F = 1 hyperfine 
ground state. The atoms are held in a single-beam optical 
dipole trap with beam waist 56 |xm, which sets the min- 
imum distance at which the field can be measured. The 
atom cloud itself has 20 pm lateral dimension, defining 
the transverse resolution. The atoms are probed with 
|xs duration pulses of linearly polarized light at 10 ps in- 
tervals, red detuned by 1.5 GHz from resonance with the 
F = 1 — > F' = transition of the D 2 line. Each pulse 
contains on average 10 7 photons. After passing through 
the atoms, the light pulses are detected by a shot-noise- 
limited balanced polarimeterP^ The experimental geom- 
etry is illustrated in Fig.^a). 

The detuning and photon number are chosen so that 
both probe-induced decoherence due to spontaneous 
emission and the perturbation due to tensorial light shifts 
are negligible during the measurement cycled This al- 
lows us to use the simple model described by Equa- 
tions ([!]) and ([2]) to fit the data. We note that the mea- 
surement sensitivity could be increased by using more 
photons and/or reducing th e det uning, at the cost of 
more elaborate data analysis! 13 * 19 ^ 

The initial atomic spin state is prepared via optical 
pumping with a single 5 ps duration circularly polarized 
pulse on resonance with the F = 1 —> F' = 1 transi- 
tion of the D2 line and propagating either along the trap 
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FIG. 1. (Color online) (a) Experimental geometry. PD: 
photodiode; L: lens; BS: beam-splitter; PBS: polarizing beam- 
splitter; WP: waveplate. (b) Schematic of the experimental 
sequence: the atoms are initially polarized along F z via op- 
tical pumping and then probed with a sequence of 1 its long 
pulses of light at 10 ps intervals. The resulting FID signal is 
recorded over 500 us. We then immediately re-polarize the 
atoms along F y with an orthogonal optical pumping beam 
and record a second FID signal. A single shot is thus ac- 
quired in 1 ms. We fit the two measurements together to a 
simple model to extract the magnetic field components Bi. 
Typical data from a single composite FID measurement are 
shown in (c) for an initially F z -polarized state and (d) for an 
initially F^-polarized state. 
For these data the field components extracted from the fit 
were B = (941, 310, 511) nT and the coherence time 
T c = 1-3 ms. 



axis, i.e. the z-axis, to prepare an F^-polarized state, or 
along the y-axis, to prepare an Fy-polarized state. Dur- 
ing the optical pumping, the atoms are uniformly illumi- 
nated with randomly polarized light on resonance with 
the F = 2 — > F' = 2 transition of the D 2 line to prevent 
atoms accumulating the F = 2 hyperfine state. A single 
composite FID measurement consists of first preparing 
an F z -polarized state and measuring the FID signal over 
500 ps, then immediately preparing an i^-polarized state 
and again making a FID measurement. A single shot is 
thus acquired in 1 ms. 

To illustrate the technique, we first record the labora- 
tory magnetic noise at the trap, shown in Figure [2] Laser 
cooled atoms are first loaded into the dipole trap during 
~ 2 s via a two-stage magneto-optical trap (MOT). A 
small field is applied with three pairs of Helmholtz coils, 
and the experiment is triggered on the 50 Hz signal of the 
laboratory mains line. The field is sampled at a sequence 
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FIG. 2. (Color online) Recording of low frequency labora- 
tory magnetic field noise, acquired with one sample point per 
cycle and varying phase with respect to the 50 Hz signal of the 
laboratory mains line, as in equivalent-time sampling with a 
digital oscilloscope. We plot the field components estimated 
from the three measurement rrij at each point p 4 . Inset: Three 
consecutive measurements are made at each point. Here we 
plot the B x component estimated from the three measure- 
ments. Error bars are ±la statistical errors. 



of points pi at 20 ms intervals, with a variable wait time 
T before the first point. At each pi we make three con- 
secutive composite FID measurements rrij, as shown in 
the inset of Figure [2] The entire sequence is repeated 300 
times to collect statistics. 

The results show good predictability of the field from 
one cycle to the next, with a typical statistical uncer- 
tainty of <TBi — 40 nT for each field component. Note 
that the experiment has no magnetic shielding, so that 
the observed variance is dominated by magnetic field 
noise from the laboratory environment. We observe 
T2 ~ 1.5 ms, which sets a limit on the coherence time 
and is important for design of future experiments. The 
FID signals give information about one gradient compo- 
nent, the one along the average field. With three FIDs, 
with applied bias fields along different directions, we can 
obtain dTS/dz, all the gradient components affecting the 
experiment. 

We are interested in our ability to predict (or retro- 
diet) the magnetic field at a moment shortly after (or 
before) the magnetic field measurement. This ability de- 
termines the precision of corrections for, or control of, 
the field seen by the atoms. To quantify this precision, 
we analyze the conditional uncertainty between consec- 
utive measurements mi, m% at T = 29, 30 ms, shown in 
Figure [51 inset. Typical experimental data are shown in 
Figure pi For a single parameter, the conditional vari- 
ance is var(y|x) = var(y — ax), where the correlation pa- 
rameter a = cov(x, y)/v&r{x) minimizes the conditional 
variance. This is equivalent to minimizing the residuals 
of a linear regression y = ax + b, and is illustrated for a 
single parameter in Figure p[b)-(e). For the data shown, 




0.2 0.4 0.6 0.8 1.0 

Sy,m, G"T) 

FIG. 3. (Color online) (a) Correlation plot of the field com- 
ponents Bi estimated from consecutive measurements mi and 
77i2 of the field. Marginal distributions from the B x compo- 
nent (b) are shown in histograms (c) and (d). Each has a 
standard deviation of ovb^ — 40 nT. Solid lines are gaussians 
with the indicated mean and standard deviation, suggesting 
that the field noise is approximately normally distributed. 
The conditional uncertainty of the measurement 7712 given the 
outcome of measurement mi is equivalent to the dispersion 
of the residuals of a linear regression of 7712 on mi, shown in 
inset (e), and has a standard deviation cr m2 | mi = 9 nT. 

the conditional uncertainty is cr m2 i mi = 9 nT. 

This analysis is readily extended to multivariate data. 
If x and y are vectors of parameters, with covariance 
matrices T x = cov(xi,Xj), and T xy = cov(xi,yj), then 
the conditional covariance matrix is given by 

Ty|x = Ty — Ty r x y . (7) 

The matrix of coefficients A = T~ F x y minimizes the 
mean squared error of the linear regression y = Ax+BP^ 
For the data shown in Figure [3j the covariance matrix 
for the first measurement is 

/ 1.60 0.21 -0.14 \ 
T x = 0.21 0.42 0.00 x 10~ 15 T 2 . (8) 
\ -0.14 0.00 0.24 J 

The corresponding conditional covariance matrix is 

/ 5.6 -8.0 -6.0 \ 
T y | x = -8.0 36.0 1.5 x 10~ 17 T 2 . (9) 
\ -6.0 1.5 12.0 J 

This shows strong correlations among the different B 
components, and it is interesting to diagonalize r y | x to 
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find uncertainties (5B 1 ,SB 2 ,SB 3 ) = (19.5, 12.0, 3.0) nT 
along the directions ei = (—0.26,0.96,0.11), e 2 = 
(-0.41,-0.22,0.89) and e 3 = (0.87,0.18,0.45), respec- 
tively. We note that e3 is nearly the field direction, indi- 
cating good predictability for the magnitude of the field. 
We observe similar results if we analyze the correlation 
between two measurements Pi at the same phase on dif- 
ferent cycles of the 50 Hz mains line. It should be noted 
that these results include readout noise, which we now 
compute. 

Faraday rotation measurement at or near the shot- 
noise limit has been demonstrated with a variety of 
cold atom systems, including releas ed MO TsP^H optical 
lattices,^ and optical dipole trapsl^^H Our experi ment 
is shot-noise limited by 10 dB at 10 7 photons/pulsePSEIl 
We compute the shot-noise-limited sensitivity using 
Fisher Information (FI) theory. 2 ®-' For a normally- 
distributed random variable 9 with fixed variance a 2 
and mean 0(x), where x is a vector of parameters, the 
FI matrix is = <r~ 2 [di8]dj9, where d a represents 
d/dx a . This directly gives the covariance matrix for x as 
T x = I . Due to shot noise, the measured rotation an- 
gles are normally distributed with a 2 . — 1/N P and means 
0i,2 from Eqs. 0, Also, the FI is additive over 

independent measurements, so the FI matrix from FID 
is 

x a = N^^MtkWMtk)- (io) 
fc,i 

where x = [B x , B y , B z , T 2 , F z (0) , F y (0)] and {t k } are the 
measurement times. 

Considering 7 = 2tt x 7 GHz/T for the ground states 
of 87 Rb, and typical values from the data of Figure [3] 
(B x ,B y ,B z ) = (910,285,540) nT, T 2 = 1.48 ms, F z (0) = 
F y (0) = 10 5 spins, G = 0.89 x 10~ 7 rad/spin and N p = 
10 7 photons, we find the covariance matrix (B portion 
only) 

/ 1.30 -2.43 -1.00 \ 
r SN = -2.43 11.87 -1.61 x 10~ 17 T 2 . (11) 
\ -1.00 -1.61 2.57 J 

If we diagonalize Tsn we find uncertainties 
{5B 1 ,5B 2 ,5B 3 ) = (11.2,5.6,0.6) nT, along the di- 
rections ei = (0.2, -0.97, 0.14), e 2 = (0.5,-0.02,-0.86) 
and e 3 = (0.84,0.24,0.48), respectively. 

We can now correct the measured field noise ofEq. f9} 
for the measurement noise, to find the field noise Tb = 
r y | x — Tsn of 

/ 4.30 -5.50 -5.00 \ 
T B = -5.50 24.00 3.10 x 10~ 17 T 2 (12) 
\ -5.00 3.10 9.50 / 

or SBi ~ 10 nT integrated over the kHz bandwidth of the 
measurement. 

The FI analysis also reveals that SF z (0) and 8F y (Q), 
the noises in the atomic state preparation, are only very 



weakly coupled into the estimates of B and T 2 , making 
the measurement insensitive to, e.g., atom number fluc- 
tuations and variation in the optical pumping efficiency. 

We have demonstrated a cold-atom magnetometer 
with sub-nT sensitivity, 20 \im transverse spatial reso- 
lution and 1 ms temporal resolution. The instrument 
gives simultaneous information about the three field com- 
ponents plus one gradient component and requires no 
additional applied fields, making it very attractive for 
non-disturbing field monitoring and control. We note 
that sensitivity can be improved by increasing the num- 
ber of atoms (in our system a five-fold improvement to 
10 6 atoms is readily achievable*^) and/or the number of 
photons, although tensor light shift s should be taken into 
account for larger photon numbers! 13 * 22 * 29 ^ 
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